Mermin-Ho vortex in ferromagnetic spinor Bose-Einstein condensates 
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The Mermin-Ho and Anderson- Toulouse coreless non-singular vortices are demonstrated to be 
thermodynamically stable in ferromagnetic spinor Bose-Einstein condensates with the hyperfine 
state F — 1. The phase diagram is established in a plane of the rotation drive vs the total mag- 
netization by comparing the energies for other competing non-axis-symmetric or singular vortices. 
Their stability is also checked by evaluating collective modes. 

PACS numbers: 03.75.Fi, 67.57.Fg, 05.30.Jp 
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Topological structure plays an important and decisive 
role in various research fields, ranging from condensed 
matter physics to high energy physics. They provide a 
common framework to connect diverse fields, enhancing 
mutual understanding jjj. 

Recent advance of experimental techniques on Bose- 
Einstein condensation (BEC) prompts us to closely and 
seriously look into theoretical possibilities which were 
mere imagination for theoreticians in this field. This 
is particularly true for spinor BEC where all hyperfine 
states of an atom Bose-condensed simultaneously, keep- 
ing these "spin" states degenerate and active. Recently, 
Barrett at al ^ have succeeded in cooling 87 Rb with the 
hyperfine state F — 1 by all optical methods without 
resorting to a usual magnetic trap in which the inter- 
nal degrees of freedom is frozen. Since the spin interac- 
tion of the 87 Rb atomic system is ferromagnetic, based 
on the refined calculation of the atomic interaction pa- 
rameters by Klausen at alB, we now obtain concrete 
examples of the three component spinor BEC (F = 1, 
mp = 1,0,-1) for both antiferromagnetic ( 23 Na)Q and 
ferromagnetic interaction cases. In the present spinor 
BEC the degenerate internal degrees of freedom play 
an essential role to determine the fundamental physi- 
cal properties. There is a rich variety of topological de- 
fect structures, which are already predicted in the ear- 



lier studies 
by others 




on the spinor BEC. These are followed 
§ § @ § @ § who exam- 
ine these topological structures more closely, such as 
skyrmion, monopole, meron or axis-symmetric or non 
axis-symmetric vortices both for antiferromagnetic and 
ferromagnetic cases. 

Superfluid 3 He is analogous to the spinor BEC where 
the neutral Cooper pair possesses the orbital and spin 
degrees of freedom, thus the order parameter is a 
multi-component |r7| . Among various topological struc- 
tures, the Mermin-Ho (MH)Jl8| and Anderson- Toulouse 
(AT) pl| vortices of a coreless and non-singular ^-vector 
texture are proposed in 3 He-A phase. These are an ex- 
tremely interesting object to study if they exist. The MH 
vortex is expected to spontaneously appear in a cylindri- 
cal vessel without any external rotation as an equilib- 



rium state because the rigid wall boundary imposes the 
Z-vector perpendicular to the vessel wall. The MH vor- 
tex is stable also under slow rotation because of their 
non-singular coreless structures 

A similar topological structure, called skyrmion in 
general is proposed in the spinor BEC. Khawaja and 
Stoofjio) study a skyrmion in the ferromagnetic BEC 
with F = 1 and conclude that the skyrmion is not a 
thermodynamically stable object without rotation. Once 
created, its radius shrinks to vanish in spite of ingenious 
proposals JTl], [T^ , |l3| , |l4| as to how to created it and to 
detect it. 

Here we demonstrate that in the ferromagnetic spinor 
BEC with F = 1 trapped in a harmonic potential the 
Mermin-Ho vortex and Anderson- Toulouse vortex are 
thermodynamically stable for a certain region of the ex- 
ternal rotation frequency and the magnetization of the 
system. These vortices are stabler than the singular 
vortices and other possible non-axis symmetric vortices. 
This comes about because as a general tendency the fer- 
romagnetic interaction leads to spatial phase separation 
of the three components (</>i, 4>q, 4>-i). Then by assigning 
different winding numbers (0,1,2) to these components 
in this order, they can be arranged so as to be effec- 
tively phase-separated in the radial direction, resulting 
in a concentric layered structure under a given magne- 
tization. The central region is occupied by 4>i, (f>o is in 
the intermediate region, and 0-i is in the outer layer. 
This ingenious non-singular coreless vortex acquires the 
angular momentum to effectively lower the total energy 
under rotation. 

We start with the standard Hamiltonian by Ohmi and 
Machida§, and Ho§: H rot = H - Jdrft ■ £\ *](r x 
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and i,j,k,l = (0,±1) corresponding to the above three 
species. The chemical potentials for the three com- 
ponents; jjii (i = 0, ±1) satisfy /i X + = 2/x . 
The scalar field V(r) = ^m(2nv r ) 2 (x 2 + y 2 ) is the 
external confinement potential such as an optical po- 
tential. The scattering lengths clq and ct2 charac- 
terize collisions between atoms with the total spin 
and 2 respectively. As mentioned, the recent refined 
estimate^] for 87 Rb concludes it ferromagnetic (do > 
02) and g s /g n = —0.01 ~ —0.005. The external 
rotation frequency Q is normalized by the harmonic 
confining frequency. It readily leads to the Gross- 
Pitaevskii (GP) equation extended to the three compo- 
nent order parameters: [{/i(r) — /x, + g n J2k \4>k\ 2 } % + 
SsE«Efei{(^)ii( f <.)iti^!}Wj = °- These coupled 
equations for the j-th condensate wave function <j)j — 
(i&j) (j — 0, ±1) are used to calculate various properties 
of a vortex in the following. 

It is important to realize that in a real experimental 
situation the total number N — J dr£j|^j| 2 and the to- 
tal magnetization M = j drSj j |0j | 2 are fixed. In our 
case these quantities are fixed by adjusting the chemi- 
cal potential (/io) and the fictitious magnetic field (// = 
Mi — Mo). The extended GP equation is solved numer- 
ically with two different methods for a two-dimensional 
disk: One is not to assume axis-symmetry. This calcu- 
lation is backed up by the computation where the axis- 
symmetry is assumed. The actual calculations are done 
by discretizing two-dimensional space into 51x51 mesh. 
The vortex configuration is characterized by the combi- 
nation of the winding number of each component. We 
denote it as (uii, Wo, W-\) where integers W\,wq, W-i are 
the winding numbers of <j)\,<j)o,<j)-\ respectively where w 
means the phase change by 2irw when the wave function 
goes around the phase singularity. 

We have performed extensive search to find a sta- 
ble vortex, starting with possible vortex configurations, 
covering a wide range of the ferromagnetic interaction 
strength g s /g n = ~ —0.02 and examining various 
axis-symmetric and non-axis-symmetric vortices (see the 
classification of possible vortices in the axis-symmetric 
case[p5[ [l6|]). We use the following parameters: the mass 
of a °^Rb atom m=1.44x 10 _25 kg, the trapping frequency 
j/ r =200Hz, and the area density n z =2.0x 10 3 /^m. The 
results displayed here are g s /g n — —0.02. 

The Mermin-Ho (MH) vortex 

is described by (</>i, </>o, </>-i) = 
^/n(r) {cos 2 f , V2e^ sin § cos § , e 2 ^ sin 2 § ) where 
the bending angle (3{r) is < (3(r) < n, (f> is the 
polar angle in polar coordinates. The spin direc- 
tion is denoted by the so-called /-vector Q given by 
l(r) — z cos [3 + sin /3 (cos 0a; + sin^y) where (5 varies 
from ,3(0) = to 0(R) = § for MH and to 0(R) = it 
for AT (R is the outer boundary of the cloud). Thus 
the spin moment is flared out to the radial direction 
and at the circumference it points outward for MH and 




FIG. 1: Density profile of the condensates (a) and the density 
map of / z -vector (b) for the (0, 1,2) vortex at fl = 0.35 and 
M/N=0.21. The bold line is the total density £\ \<f>j\ 2 and 

the thin lines show the density of each component \<f)j\ 2 - 

downwards for AT (see for schematic /-vector structure 
Fig. 18 in Ref. |l7| ). These vortices have the winding 
number combination (0,1,2) in our notation. 

We show the calculated MH vortex in Fig. 1 where the 
axis-symmetric density profiles for each component and 
the density map of l z are displayed. It is seen that the 
central region of the harmonic trap is occupied by <pi 
with zero winding number w± — 0. The 0o component 
which has a singularity with wq = 1 at r = is pushed 
outward while the </>_i component occupies the circum- 
ference region because of W2 — 2. Note that in the small 
r region the singular vortex with w behaves like oc r w . 
The resulting total density is coreless and non-singular. 
It has a smooth density variation described by a Gaussisn 
form except for the outermost region. 

In Fig. 2 we show the spatial dependence of the l z - 
component along the radial direction, namely, the spatial 
dependence of the bending angle (3(r) for the MH vortex. 
As the magnetization M decreases, the local magneti- 
zation changes from positive to negative values through 
zero. It means that the Z-vector in the vortex flares out 
radially to orient almost horizontally (3(r = R) = ^ for 
M/N ~ 0.5 and to point down for (3(r = R) = tt for 
M/N ~ 0. The former (latter) corresponds literally to 
the Mermin-Ho (Anderson- Toulouse) vortex. This is sim- 
ply because as M decreases, the spin-down component 
0_i with w = 2 increases in the outer region. Thus we 
can control these MH and AT vortices by merely chang- 
ing the total magnetization. We notice that in superfluid 
3 He-A phase the stability of the MH and AT vortices is 
guaranteed by the boundary condition at the wall of a 
vessel where the Z-vector is constraint to be parallel to 
the wall surface fTl. Here the situation is completely dif- 
ferent; we impose no constraint on the /-vector direction. 
These vortex configurations are created naturally under 
the condition of a given total number and magnetization, 
both of which are well controlled in a harmonic trap ex- 
periment. 

This difference in the two systems is rather remark- 
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FIG. 2: Spatial dependence of the Z z -component along the 
radial direction at M/N=0, 0.46, and fl = 0.37. The dashed 
line shows cos/3(r) with the bending angle P(r) = nr/R (R = 
2.85fim) for the AT vortex. 



able and interesting. In our case MH or AT vortex comes 
about purely due to the spin interaction, which is writ- 
ten as oc \g s \\2cf)i(r)(f)^i(r) — (j>o(r)\ 2 . This term favors 
the mutual phase segregation jl5j |lq |. If 4>\{r) ^ and 
4>-i (r) ^ 0, 4>o(r) comes in, explaining the concentric lay- 
ered structure. In fact, in the antiferromagnetic case the 
similar calculation shows that the MH and AT vortices 
are never stabilized under a similar condition Jl^, gjj . 

We have done extensive search for determining the re- 
gion for the (0,1,2) vortex of the MH and AT in the 
plane of f2 vs M/N . We check our computation by start- 
ing with a certain initial vortex configuration, allowing 
non-axis-symmetric vortex. We end up with the axis- 
symmetric MH and AT vortices where these are stable. 
The stability for these vortices is also examined from the 
two aspects: one is the global stability and the other is lo- 
cal stability in the energy landscape. The global stability 
means comparing the relative energy of various vortices 
and vortex free state to select out the lowest energy state. 
The local stability is discussed shortly. 

The resulting phase diagram is displayed in Fig. 3 
where a large area is occupied by the (0,1,2) vortex, in- 
cluding MH and AT. Near M ~ the non-axis-symmetric 
(1, 1, 1) vortex and near M/N ~ 1, (1,0, —1) vortex are 
stabilized. We find a large empty region in the interme- 
diate M/N region where no vortex and vortex-free state 
are stabilized at all because the phase separation in the 
ferromagnetic case prevents forming a uniform mixture 
of the three components even when the circulation is ab- 
sent for the vortex-free state. This is in contrast with the 
antiferromagnetic case in which everywhere is occupied 
by a stable phase. 

It is noted that the focused frequency region = 0.38 
corresponds to the single vortex region in the scalar BEC 
case[^2| beyond which multiple vortex state must be se- 
riously considered, thus the present single vortex consid- 
eration is justified in the lower frequencies. 

We examine various competing vortices, including a 
singular vs non-singular and axis-symmetric vs non- 
axis-symmetric vortices. In particular, since for axis- 
symmetric vortices the winding number combination of 
the three components (wi,wq, w-i) is restricted to 2wq = 
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FIG. 3: Phase diagram in the ferromagnetic state. The 
lowest eigenvalues along qe at 57=0.3 are displayed in the inset. 
Triangles: M/N=0.2 and filled triangles: 0.75. 



W\ + we exhaust all the possible vortices with 

Wi, Wo, W-i smaller than unity[jl6|. Note that MH and 
AT belong to this axis-symmetry category. 

The (1, 1, 1) vortex shown in Fig. 4 is non-singular and 
non-axis-symmetric which is stable in low M/N. This 
(1,1,1) vortex is advantageous because (A) it does not 
contain the higher winding number. It generally leads to 
collapse fewer winding number vortex, such as 2 — ► 1 + 1 
vortices. (B) The overall condensation energy is gained 
by placing their singularities off the trap center where 
the potential energy is minimum. (C) As Q increases, 
the two separate singularities of <j)\ and </>_i adjust their 
mutual distance from the center so as to maximize the 
angular momentum L in order to save the energy — Q ■ L. 
In this sense this configuration is flexible against varying 
fl. These reasons explain why this vortex survives along 
the ft-axis near M/N ~ in Fig.3. 

In contrast, the (0,1,2) vortex contains the higher 
winding number 2 for which makes it less advanta- 
geous against the (1, 1, 1) vortex. However, upon varying 
M , this vortex is quite flexible by adjusting the particle 
numbers for the three components; As M increases, the 
number of <f>\ grows smoothly relative to the rests. The 
non-winding <pi component works as a "pinning" center 
for the remaining 0o and <p-i- In particular, 0_i with 
W-i = 2 is stabilized by the presence of <j>\. Therefore 
as M increases or <pi grows, the MH and AT vortices 
become more and more stable. This explains that the 
critical frequency, or the lower phase boundary of MH in 
Fig.3 becomes less as M increases. As for the (1, 1, 1) 
vortex, in comparison, is not flexible enough for varying 
M because the symmetric arrangement of the two sin- 
gularities in <f>i and <p_i around the center (see Fig. 4) 
becomes unbalanced as <f>\ component grows while 
shrinks. Thus it is confined in a narrow thin region near 
M/N - (see Fig.3). 

The (1,0,-1) vortex, (not shown here, see Fig.3 in 
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FIG. 4: Density profile of the condensates for the non-axis- 
symmetric (1,1,1) vortex in Q = 0.35 and M/N=0.0. The 
bold line is the total density £ I0j| 2 an d the thin lines show 
the density of each component \<j)j\ 2 - 



Ref. JT5| ) is stabilized in larger M/N and higher f2 region 
because in this vortex the dominant component </>i has 
the winding number 1, which can effectively lower the 
rotation energy. 

Having finished the "global" stability of the MH and 
AT vortices, we now turn to the "local" stability, that 
is, the stability against a small perturbation. This is 
done by solving the extended Bogoliubov equations 
to the three components under the axis-symmetric 
situation HI): J2j{ A ij u q( r J) ~ B ij v q( r J)} = £ q u q (v,i) 



and £j{£tf« 9 (r.i) 



ij = h(r)Sij — 

.EaEfcjK^Jy^aJw^l + (F a )u(F a ) kj (tlfr)}, 
9n^i4>j + 9sT, a J2ki [(F a )ik(j)k(F a )ji(j)i} , u q (r,i) 



Bij - 



and v q (r,i) are the g-th eigenfunctions with the spin i 
and e q corresponds to the q-th. eigenvalue. Since this 
gives the excitation spectrum of the system, the negative 
energy relative to the condensation energy at zero im- 
plies the local intrinsic instability of the putative vortex 
in the energy landscape. We have performed the exten- 
sive computation to check this local stability for MH 
vortex and other axis-symmetric vortices against a small 
perturbation. As expected, the lower phase boundary of 
the MH in Fig. 3 coincides almost completely with the 
local stability region, below which the lowest excitation 
mode with the angular momentum qg = — 1 for the 
positive external rotation (f2 > 0) becomes negative (see 
the inset in Fig. 3). The upper boundary of the (0, 1,2) 
vortex in Fig. 3 indicates that the lowest excitation mode 
with mainly qg = +1 becomes negative (see the inset 
in Fig. 3). Thus the global stability mentioned roughly 
corresponds to the present local stability. Therefore, it 
is concluded that the MH and AT vortices are stable 
and robust in the ferromagnetic spinor BEC. 

It is easy to calculate the total angular momentum 
L in MH vortex which is given by Ljfi = Nq + 27V_i 
(Ni is the particle number of the i-component). Since 
M = Ni—N^i and N = Ni+No+N-i, the total angular 
momentum per particle is found to be ^ = 1 — jq-. This 
simple formula means that at M = 0, <fio with wo — 1 
carries all the angular momentum as in the usual scalar 



vortex, and at M = A, L = because (f>i has no winding. 
At 4£ = -g, jlfj is half, implying that in the MH vortex 
the angular momentum per particle is exactly h/2. 

In summary, we have shown that the Mermin-Ho and 
Anderson- Toulouse vortices are thcrmodynamically sta- 
ble under a certain rotation drive and the total magne- 
tization of the system, both of which are experimentally 
well-controlled parameters. We have examined their sta- 
bility by two methods, local and global ones in the energy 
landscape, coinciding with the earlier conclusion that 
these are unstable under no rotation drive |^J. These in- 
triguing objects might be detected by various ways, such 
as optically by utilizing the Faraday rotation which im- 
ages the spatial magnetic pattern of these vortices. A fa- 
vorable experiment is to use an oblate ellipsoidal shaped 
system to mimic our disk calculation. An elongated cigar 
system may not be appropriate, which induces the phase 
separation along the long axis. 

We thank T. Ohmi and T. Isoshima for valuable dis- 
cussions. 
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